Crystal-field ground state of the non-centrosymmetric superconductor CePt 3 Si: 
a combined polarized soft X-ray absorption and polarized neutron study 
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We determined the crystal-field split Hund's rule ground state of the non-centrosymmetric heavy 
fermion superconductor CePtaSi with polarization dependent soft X-ray absorption spectroscopy 
(XAS) and polarized neutron scattering. We are also able to give the sequence of the crystal-field 
states from the temperature evolution of the linear dichroic signal in the XAS. The quantitative 
analysis of the XAS temperature dependence together with the neutron transition energies complete 
the identification of the crystal-field level scheme. 

PACS numbers: 71.27.+a, 75.10.Dg, 78.70.Dm, 78.70Nx 
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I. INTRODUCTION 

The discovery of superconductivity in the antifer- 
romagnetic phases of the non-centrosymmetric heavy 
fermion compounds CePtaSi, 1 CeIrSi 3) 2 and CeRhSi 3 3 
stimulated a flurry of theoretical and experimental 
activitieSM^jl with the main interest to understand the 
pairing state in the absence of inversion symmetry. It 
seems that the superconductivity is due to a mixture of 
singlet and triplet Cooper pairing and that the super- 
conducting state is magnetically driveni&£ The unusual 
properties evolve from the hybridized 4/ electrons so that 
the question arises which influence the highly anisotropic 
crystal-field ground state may have on the hybridization 
between 4/ and conduction electrons. Attempts were 
made to investigate this issue for the heavy fermion com- 
pound CeRu 2 Si2, 10 and the semiconducting Kondo mate- 
rials CeNiSn ) 11 ' 12 and CeRhAsi 3 .. More recently the as- 
pect of g-dependent hybridization was also studied in the 
CeMIns compounds . 14 ' 15 ' 16 ! 17 However, reliable experi- 
mental input for such theoretical investigation is scarce 
due to experimental uncertainties in determining the 4/ 
ground state wave functions. 

The wave functions and crystal-field levels are usu- 
ally determined from a combined analysis of single crys- 
tal susceptibility and polycrystalline neutron scattering 
measurements. However, the method often fails because 
of thermal or powder averaging in either technique, and 
additionally, neutron scattering suffers from phonon con- 
tributions and broadened crystal-field levels due to hy- 
bridization effects*^ Moreover, in the presence of mag- 
netic order, anisotropic molecular field parameters have 
to be introduced in order to fit the static susceptibility, 
ending up with too many free parameters for a unique 
description. The methods to apply are then inelastic or 
elastic polarized neutron scattering on single crystals, the 
first technique requiring large single crystals, the second 



being very sensitive to absorption corrections when de- 
termining the 4/ magnetic form factor. 

We have shown for the case of the heavy fermion mate- 
rial CePd2Si2 that soft X-ray absorption at the Ce M^s 
edges can be used complementary to neutron scattering 
in order to determine the ground state wave function. 19 
XAS is highly sensitive to the initial state and via its po- 
larization dependence direct information about the \J Z ) 
admixtures of the ground state wave function can be ob- 
tained. Sensitivity for higher lying crystal-field states 
is achieved by thermally populating those states. It 
was shown in Ref. that the capability of XAS to re- 
solve closely lying crystal-field states is not determined 
by the energy resolution of the experimental set-up, but 
rather by the temperature. However, the exact energy 
of crystal-field excitations is most accurately determined 
with inelastic polarized neutron scattering on single crys- 
tals where phonon and magnetic scattering can be sepa- 
rated. Hence it is ideal to combine both techniques. 

CePtaSi orders antiferromagnetically at Tjy=2.2 K and 
enters a superconducting state at T sc «0.75 The co- 
existence of both states has been confirmed by several 
experimental techniques ! 20 ' 21 The low-energy magnetic 
excitations were recently measured on a single crystal us- 
ing inelastic neutron scattering. 22 CcPtaSi crystallizes in 
the non-centrosymmetric tetragonal space group P4mm 
with lattice constants a = 4.072 A and c = 5.442 AM 
Inversion symmetry is broken by the absence of a mirror 
plane perpendicular to the c-axis. However, the lack of 
inversion symmetry is irrelevant for the crystal-field de- 
scription since the 4/ wave functions are even functions. 
Therefore, as in other centrosymmetric compounds, the 
Hund's rule ground state of Ce 3+ (J = 5/2, point group 
C^ v £ Dm)— splits under the influence of a tetragonal 
crystal- field, into three Kramer's doublets, which can be 
represented in the basis of |J 2 ). The eigenfunctions of 
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the three Kramer's doublets can be written as 
|2) = | ±1/2) 

|1> = /3|±5/2}- a | T 3/2) (1) 
|0) = a|±5/2}+/3| T 3/2) 

with a 2 + f3 2 — 1 and an arbitrarily chosen phase of the 
mixed states. 

Two groups have published inelastic neutron scattering 
data of CePtaSi powder samples in order to determine 
the crystal-field scheme. Bauer et al3£ found a broad 
distribution of magnetic intensity between 10 meV and 
30 meV which was fitted with two non-resolved crystal- 
field excitations centered at 13 meV and 20 meV. The 
crystal-field parameters of Bauer et al. yield |0) = 
0.59| ± 5/2} + 0.81| T 3/2) for the ground state wave 
function and |2) = | ±1/2) as the highest excited state at 
20 meV. Metoki et alj2& found also a mixed ground state, 
but interpreted their neutron data with a very different 
ground state mixing of |0) = 0.93| ± 5/2} + 0.37| T 3/2). 
According to their data the | ± 1/2) state is 1 meV above 
the ground state and the second mixed state is at 24 meV. 

In this paper, we determine the ground state wave 
function from the low temperature XAS data, i.e. the 
value of the mixing factor a, which gives the actual spa- 
tial distribution of the crystal-field ground state. We 
will further show that the temperature dependence of the 
XAS data determines the sequence of states. A quantita- 
tive modeling of the temperature dependence can even be 
done provided that the crystal-field energies are known. 
The latter are determined from the polarized inelastic 
neutron data presented here. Thus a complete picture of 
the crystal-field level scheme has been achieved. 

II. EXPERIMENTAL 

The CePtsSi single crystal used for this study was the 
same as in Ref. [22|. A thin slice was cut from this crys- 
tal and used for the XAS measurements. We recorded 
all spectra using the total electron yield method in a 
chamber with a pressure of 4 x 10~ n mbar at the UE- 
46 PGM-1 undulator beam line of BESSY II. The total 
electron yield signal was normalized to the incoming pho- 
ton flux Iq as measured at the refocusing mirror. Clean 
sample surfaces were obtained by cleaving the samples 
in situ. The energy resolution at the cerium M4 5 edges 
(hu « 875 - 910 eV) was set at 0.15 eV. Further/the un- 
dulator together with a normal incident measurement ge- 
ometry allow for a change of polarization without chang- 
ing the probed spot on the sample surface. This guaran- 
tees a reliable comparison of the spectral line shapes. The 
two polarizations were E ± c and £7 1| c, c being the long 
tetragonal axis. The sample was cleaved and then cooled 
down to the lowest accessible temperature. Spectra were 
taken at 5 K, 50 K, 150 K, and 200 K while warming up. 
In order to assure reproducibility of the data the sample 
was cleaved a second time and then data were taken in 
the sequence 300 K, 50 K and 80 K. 



Inelastic neutron scattering measurements using po- 
larization analysis were performed on the thermal triple- 
axis spectrometer IN20 at the high-flux reactor of the 
Institut Laue-Langevin. Horizontally focusing Heusler 
monochromator and analyzer were used in combination 
with relaxed collimation and a fixed final energy of 14.7 
meV with an energy resolution (full width at half max- 
imum, FWHM) of 0.6 meV. Higher-order neutrons were 
removed from the scattered beam using a graphite fil- 
ter and their contribution to the monitor count rate was 
corrected for. Measurements were performed in the para- 
magnetic phase at a temperature of 5 K with the wave- 
vector transfer Q either parallel to the a or to the c axis, 
but with the same value of Q = |Q|, to eliminate the 
influence of the magnetic form factor to the measured in- 
tensities. The use of polarization analysis in combination 
with a single crystal has two advantages for crystal-field 
determinations compared with standard methods of un- 
polarizcd neutrons on a polycrystalline material: (i) Scat- 
tering from phonons can easily be subtracted; (ii) Sepa- 
ration of the j = x, y, z components of the angular mo- 
mentum operator Mj to the transition matrix elements 
\M±\ 2 = \ {f\Mj-\i)\ 2 between the initial state i and the 
final state /, where ± reflects the fact that only compo- 
nents perpendicular to Q are observed in neutron scat- 
tering experiments. 

We measured the three polarization components 7 = 
(£,77,0 of the spin-flip cross-section <r 7 , where £ is par- 
allel to Q, 77 is perpendicular to Q in the horizon- 
tal scattering plane, and C is vertical. The pure mag- 
netic scattering is given by the difference of the spin- 
flip cross-sections || Q and 1 Q or, more specifically for 
the present set-up, the different transition matrix ele- 
ments of the Mj operator are given by \M Z \ 2 = — <r v 
for Q || a and IM^ 2 = \M y \ 2 = - cr c for Q || a or 
\M X \ 2 = \M y \ 2 =(7 6 - a n = (T C - cr c for Q || c. The three 
independent determinations of l-M^I 2 are identical within 
statistical uncertainty. The measured non-spin-flip cross- 
section shows no strong intensity that could leak through 
to the spin-flip channels via the finite flipping ratio of 20- 
23, which corresponds to a neutron polarization efficiency 
of 95% at both polarizer and analyzer. 



III. RESULTS AND DISCUSSION 

For the C^v point symmetry the crystal-field Hamil- 
tonian H CF = B? 2 0% + B%C>1 + B\0\ describes the 
crystal field potential when the three Stevens parame- 
ters B®, B®, and B\ are determined. Alternatively the 
crystal-field potential can be expressed in terms of the 
transition energies E\ and E2 between the crystal-field 
states within the multiplct and the mixing parameter 
a. We performed ionic full multiplct calculations us- 
ing the XTLS 8.3 program 2 ^ to calculate the XAS spec- 
tra. All atomic parameters were given by Hartree-Fock 
values, with a reduction of about 35% for the 4/ — 4/ 
Coulomb interactions and about 20% for the 3c? — 4/ in- 
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FIG. 1: (color online) Measured and simulated linear polar- 
ized XAS spectra of CePtsSi at the IVU^ edges at T—5 K. 
The inset shows the corresponding spatial distribution. 

teractions to reproduce best the experimental isotropic 
spectra, Jj so = 2ij_ + In . This accounts for the config- 
uration interaction effects not included in the Hartree- 
Fock scheme. Once the atomic parameters are fine tuned 
to the isotropic spectra, the polarized XAS data can be 
described by the incoherent sums of the respective po- 
larization dependent spectra of the pure \J Z ) states^ as 
long as the crystal-field splitting is small with respect to 
the spin orbit splitting. The latter requirement, which 
is fulfilled here (E so m 280 meV and E CF < 20 mcV) 
assures that interference terms resulting from intermix- 
ing between the J = 5/2 and J — 7/2 multiplet can be 
neglected. 

Figure 1 shows the low temperature linear polarized 
XAS data of CePtaSi at the M^g edge and their simu- 
lation based on the full multiplet treatment as described 
above. The polarization effect is not very large but real 
as will be shown in the next paragraph from the temper- 
ature dependence of the XAS data and the reproducibil- 
ity of the 50 K spectra after the second cleave. The pure 
| ± 1/2} state cannot describe our low temperature data 
as can be seen from Fig. 2, where we have shown again 
(see Ref. H3) for reasons of convenience, the polarization 
dependent spectra for each \J Z ) state. The 5 K data are 
well described with the mixed ground state 

|0) =0.46|±5/2) + 0.89|t3/2) 

(a = 0.46 ± 0.01) without the necessity to include the 
fractional occupation of a higher lying state. The small 
size of the polarization effect can be understood when 
considering, as was shown in the paper by Hansmann 
et al.fi^ that the variation of a from to 1 leads to a 
vanishing polarization effect when a is equal to y/l/6 — 
0.41, i.e. when a is equal to the value of a cubic crystal- 
field symmetry. According to our data the mixing factor 
in CePtaSi is close to this value. 

In order to rule out that the small size of the effect 
is due to unwanted surface effects the sample was mea- 
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FIG. 2: (color online) Simulated polarizations for pure \J Z ) 
states of Ce 3+ and sum of the three polarization dependent 
\J Z ) spectra. The insets resemble the spatial distributions. 



sured at several temperatures, cleaved a second time and 
remeasured. The left panel side of Fig. 3 shows the M5 
edge of CePtaSi for several temperatures between 5 K 
and 300 K. After the first cleave data were taken at 5 K, 
50 K, 150 K and 200 K. After the second cleave the sam- 
ple was measured first at 300 K and then cooled down 
to 50 K and warmed up to 80 K. Comparing the data 
at different temperatures shows that the polarization in- 
creases from 5 K to 50 K and 80 K and then decreases 
again as temperature rises further. We conclude that this 
effect is real since the polarization at 300 K (measured 
first after second cleave) is smaller than the polarization 
at 50 K (measured second after second cleave) and 80 K 
(measured third after second cleave). Surely, if surface 
degradation with time had been present, then one would 
rather expect to see a decrease of polarization. In addi- 
tion the 50 K data were well reproduced. The black dots 
in Fig. 4 are the difference of the two intensities E _L c 
and E || c (called linear dichroism LD) at the high energy 
peak A of the M 5 edge (see Fig. 3) and they summarize 
the temperature evolution of the polarization effect. 

Polarized inelastic neutron scattering was used to de- 
termine the energy of the crystal-field levels. Figure 5 
shows energy scans for Q=(1.5,0,0) for difference cross- 
sections sensitive to M x and M z , respectively. Within 
the crystal-field scheme of Eq. (fTJ) , the M z cross-section 
(open blue symbols in Fig. 5) contributes only to the 
transition between the ground state |0) and the first ex- 
cited state |1). The data were fitted by a damped har- 
monic oscillator (DHO) 28 convoluted with the instrumen- 
tal resolution function. The first excited state |1) has 
an energy (with respect to the ground-state energy) of 
£i=15.0(4) meV and a width (FWHM) of 8.2(9) meV, 
which reflects the damping of the crystal-field excitation 
due to the hybridization with the conduction electrons. 
The M x cross-section contains both transitions, |0) — > |1) 
and |0) — > 1 2), and the corresponding data (full red sym- 
bols in Fig. 5) were modeled with two DHOs. The en- 
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FIG. 3: (color online) Left panel: temperature dependent XAS data at the M5 edge. Two middle panels: simulation of spectra 
with ground state wave function as determined from the 5 K XAS data and with crystal field energies obtained from the 
polarized neutron data. The simulations in the middle columns have different sequences of states. Right panel: simulation with 
ground state wave function and energies as given by Metoki et al.— 
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FIG. 4: (color online) Difference of the intensities _E_Lc and 
E || c (linear dichroism LD) at peak A of the M5 edge versus 
temperature. The black dots are determined from the experi- 
mental data (left panel in Fig. 3), the dark yellow line results 
from a simulation according to the left, middle panel in Fig. 3. 



ergy and width of the first transition (DHO) in the M x 
data were kept fixed to the above values while the inte- 
grated intensity was fixed to the fitted value of the M z 
cross-section after scaling with the corresponding transi- 
tion probabilities, calculated from the mixing parameter 
a of the XAS measurement. The second DHO described 
the second transition, for which the best fit gave an en- 
ergy of #2 = 18.7(5) meV and FWHM = 6.2(8) meV. The 



fitted integrated intensity of the second DHO is about 
16% smaller than expected from the calculated transi- 
tion probabilities. We consider this to be an excellent 
agreement between neutron intensities and our findings 
from XAS. 

The temperature dependence of the XAS data can now 
be modeled quantitatively using the ground state wave 
function as determined from XAS at T= 5 K and the en- 
ergy levels as determined with polarized neutrons. The 
results of these calculations are shown in the two middle 
panels of Fig. 3. The simulations differ in the order of the 
first and second excited state. Only the simulation with 
the I ± 1/2) as highest state reproduces the temperature 
dependence of the data well. The other option with the 
I ±1/2) as first excited state leads to an inversion of polar- 
ization as temperature rises: for example, the intensity 
of peak A for E _L c (blue) is simulated to become smaller 
than for E \\c (red) at 150 K, in contradiction with the 
experiment, see dark yellow circles in the left and middle 
right panel of Fig. 3. This can be understood when look- 
ing at the pure states as displayed in Fig. 2 and keeping 
in mind that the XAS spectra at finite temperatures re- 
semble the superimposed polarizations of each state. The 
top of Fig. 2 shows that the full thermal population of 
all crystal-field states leads to a vanishing polarization 
since the sum of all \J Z ) spatial distributions is spher- 
ical. From this can be concluded that a small ground 
state polarization requires the first and second excited 
state to almost compensate each other, i.e. the other 
mixed state must have a polarization almost opposite to 
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FIG. 5: (color online) Magnetic intensities as obtained from 
neutron scattering difference spectra taken at Q=(1.5,0,0) 
and T=5 K. Open blue symbols and line: experimental data 
of the |M Z | 2 component and fit to the transition |0) — > 
Full red symbols: experimental data of the |M^| 2 component, 
which includes the transition |0) — > |f) (dotted red line, all 
parameters fixed, see text) as well as the transition |0) —* |2) 
(dashed red line, all parameters free). The full red line is the 
sum of both contributions. 



the one of the | ± 1/2) (see Fig. 2). Hence mixing in the 
other mixed state as temperature rises leads at first to an 
increase of polarization until the influence of the highest 
state, here the | ± 1/2), compensates this due to ther- 
mal population. The dark yellow line in Fig. 4 displays 
the simulated temperature dependence of the polariza- 
tion. The simulation reproduces the experimental tem- 
perature dependence quite well, including the presence 
of a maximum at 80 K. 

For completeness we have simulated the XAS spectra 
using the crystal-field model suggested by Metoki et al»2£ 
Their ground state wave function produces spectra which 



do not resemble the experimental ones (see right panel 
of Fig. 3). Moreover, their model would lead to a strong 
temperature dependence between 5 K and 50 K which is 
not present in the experimental data. 



IV. CONCLUSION 

We have performed temperature dependent linear po- 
larized XAS measurements at the edge of cerium 
in CePtaSi and inelastic polarized neutron experiments. 
The combination of both techniques gives a conclusive 
picture of the crystal-field level scheme in CePtaSi. The 
low temperature XAS data yield 0.46|±5/2)+0.89|t3/2) 
as the ground state wave function and the sequence of 
excited states can clearly be determined from the tem- 
perature dependence of the XAS data. Using the crystal- 
field transition energies as determined from the polarized 
neutron scattering data (15.0 and 18.7 meV) we were 
also able to quantitatively explain the temperature de- 
pendence. The wave-vector and polarization dependence 
of the neutron scattering intensities also confirm this se- 
quence of the crystal-field levels. We note that the result 
with having the | ± 1/2) level as the highest lying excited 
state is in excellent agreement with the findings by Bauer 
et al. and excludes the suggestions by Metoki at al. 
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